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^ ' Abstract 

' We study the effect of dielectric inhomogeneities on the interaction 

^ ■ between two planparallel charged surfaces with oppositely charged mo- 

Cf^ i bile charges in between. The dielectric constant between the surfaces is 

I assumed to be different from the dielectric constant of the two semi- 

■ infinite regions bounded by the surfaces, giving rise to electrostatic 
, image interactions. We show that on the weak coupling level the im- 
age charge effects are generally small, making their mark only in the 
second order fluctuation term. However, in the strong coupling limit, 

I . the image effects are large and fundamental. They modify the inter- 

I actions between the two surfaces in an essential way. Our calculations 

Q ■ are particularly useful in the regime of parameters where computer 

O ■ simulations would be difficult and extremely time consuming due to 

^ I the complicated nature of the long range image potentials. 

; 1 Introduction 

Processes involving electrostatic interactions often play a dominant role in 
biological and soft-matter systems [1]. In aqueous environments charges on 
macromolecular surfaces such as membranes, self-assembled micelles, globu- 
lar proteins and fibrous polysaccharides tend to dissociate and affect a wealth 
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of functional, structural and dynamical properties [1]. Due to the nature of 
the most common solvent in these systems - water - the interfaces between 
solvent and macromolecules are not only characterised by charges but also 
by a relatively large differences in respective dielectric properties. In this 
contribution we will study the effect of these dielectric inhomogeneities on 
the interactions between charged macromolecular surfaces mediated by the 
dissolved mobile counterions. 

Understanding the behaviour of charged systems starts with the proper 
analysis of the counter- ion distribution around charged object. The tradi- 
tional approach to Coulomb fluids has been the mean-field Poisson-Boltzmann 
(PB) formalism applicable at weak surface charges, low counter-ion valency 
and high temperature [2] . The limitations of and corrections to this approach 
become practically important in highly- charged systems where counterion- 
mediated interactions between charged bodies start to deviate substantially 
from the mean- field accepted wisdom [3]. One of the most important re- 
cent advances in this field has been the systematization of these non-PB 
effects based on the notions of weak and strong coupling approximations as 
championed by Netz and coworkers [2]. These two approximations allow for 
an explicit and exact treatment of charged systems at two disjoint limiting 
conditions whereas the parameter space in between can only be analysed ap- 
proximately and is mostly accessible only via computer simulations. They 
also add a transparent and useful systematization to different approaches 
tending to upgrade or supplement the traditional PB way [4] . 

Both the weak and the strong coupling approximations are based on a 
functional integral or field-theoretic representation [5] of the grand canonical 
partition function of a system composed of fixed surface charges with inter- 
vening mobile counterions, and depend on the value of a single dimensionless 
coupling parameter S. If the charge of the counterions with valency q is qeo, 
the Bjerrum length Ib = el/ATieeokT is the length which measures the dis- 
tance at which two unit charges interact with thermal energy kT (in water 
at room temperature, the value is ~ 0.7 nm) and the surface charge den- 
sity of the bounding interfaces is a, then the coupling parameter is given by 
S = 27iq^l'^a/eo ~ q^aT'"^. For later use let us also introduce the standard 
Gouy-Chapman length defined as /i = eo/27rg/scr that measures the extent 
of the counter- ion layer next to a charged surface [1]. The meaning of the 
coupling parameter can now be easily understood by considering the mean 
distance between counter-ions in the layer next to a charged interface. The 
volume available per ion is Ana? /?> = qe^fi/a, so that the mean distance be- 
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tween counterions is of the order of a = /i(3S/2)^/^. Since the height of the 
counter-ion layer is of order of the Gouy-Chapman length, it follows that in 
the weak-coupling case, defined by S ^ 1, the height of the layer is much 
larger than the separation between two neighbouring counter-ions and thus 
the counterion layer behaves basically as a 3D gas. Each counter-ion in this 
case interacts with many others and the collective mean-field approach a la 
Poisson-Boltzmann is completely justified. On the other hand in the case 
of the strong coupling limit, defined by S 2> 1, the mean distance between 
counterions a is much larger than the layer height, meaning that the counter- 
ion layer behaves as a 2D gas [6]. In this case the mean- field approach breaks 
down, each counter-ion moving almost independently from the other ones 
along the direction perpendicular to the wall and the collective effects defin- 
ing a mean-field being non-existent. The two limits are thus characterised 
by a low/high valency of the counterions and/or a small/large value of the 
surface charge density. The range of validity of both limits has been explored 
thoroughly in [6]. 

Formally it can be shown that the weak coupling limit can be straightfor- 
wardly identified with the saddle-point approximation of the field theoretic 
representation of the grand canonical partition function [7], and is reduced to 
the PB theory to the lowest order, the next one being due to harmonic fluc- 
tuations around the saddle-point that lead to a generalized zero-frequency 
Lifshitz term [8]. The strong coupling approximation has no PB-like cor- 
relates [10] since it is formally equivalent to a single particle description, 
corresponding to two lowest order terms in the virial expansion of the grand 
canonical partition function. The consequences and the forms of these two 
limits of the Coulomb fluid description have been explored widely and in 
detail (for reviews, see [2,3]). 

In what follows we propose to add another aspect of the weak-strong 
coupling dichotomy but this time modified by the presence of dielectric in- 
homogeneities at the charged boundaries. These are ubiquitous in soft- and 
bio-systems since water, being the most common solvent, has a vastly dif- 
ferent static dielectric constant (~ 80) from all matter composing (bio) 
macromolecules (~ 2 — 5). Every macromolecule - water interface is thus 
a source of a huge dielectric jump. As the presence of dielectric discontinu- 
ities leads to (Kelvin) electrostatic image interactions, it is exactly those that 
we intend to include in the analysis of the counterion mediated interactions 
between charged macromolecular surfaces. In this work we will delimit our- 
selves solely to the planparallel geometry exemplifying the case of interacting 
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planar lipid membranes. Other geometries pertinent to interactions between 
e.g. cylindrical DNA molecules will follow in a separate publication. 

Below, we will start a la Netz and proceed via the weak and the strong 
coupling limits all the while consistently incorporating dielectric discontinu- 
ities at all levels of analysis. In the weak coupling limit that we address in 
order to present a complete and consistent analysis, we will show that for this 
particular choice of geometry the mean-field term in the interaction pressure 
is not affected by image interactions; they only transpire in the second-order 
fluctuation term around the mean-field. Since this term in the weak coupling 
limit is by necessity small, so is the image effect. The opposite is true for 
the strong coupling limit. Being effectively a single particle description it is 
quite susceptible to dielectric inhomogeneities and we will show that their 
effect on the interaction pressure is substantial. It can even add qualitatively 
new features to the interaction pressure that are absent when the dielectric 
discontinuities are ignored. The analytical results that we derive are particu- 
larly a propos in the strong coupling limit since in this case even simulations 
are hampered down by the long computer times needed to evaluate infinite 
sums of dielectric images in the energy function. 



Consider a system of counter-ions with valency q confined between two 
oppositely charged walls with uniform charge distribution of surface charge 
density a. Dielectric constant in the slab between the walls is e {e.g. water) 
whereas outside it is assumed to be different and is denoted by e' {e.g. hydro- 
carbons), see Fig. [H For simplicity we set the origin of the coordinate system 
at the middle , so that the walls are positioned symmetrically aX. z = —a and 
z = a. 

Counter-ions interact via Coulomb interaction potential n(r, r') that is 
composed of the direct interaction Mo(r, r') and the (electrostatic) image in- 
teraction Uim(r, r'), being a consequence of the dielectric constant disconti- 
nuities at the boundaries of the system, thus 



The 2D Fourier transform in lateral coordinates x, y of the two interaction 



2 The model 




(1) 



4 



counter-ions 



■ 



-a 







a 



Figure 1: Counter- ions in a medium with an inhomogeneous dielectric con- 
stant confined between two charged walls separated by 2a. The space in 
between the surfaces has the dielectric constant e, while the half-spaces out- 
side the charged walls have a dielectric constant e'. The fixed charge density 
of this system is given by a6{a — z) + a6{a + z), where a is the surface charge 
density assumed to be equal at both bounding surfaces and being of the sign 
opposite to the sign of the counterions. 

potentials is (see e.g. [9]) 



Uo{Q,z; z') 
Uim{Q,z; z') 



-Q\z-z>\ 



2eeoQ 

1 chQ(z + z') + Ae-2««chQ(z-^') 



eeoQ 



^-lg2Qa _ ^ g-2Qa 



(2) 



where Q is the magnitude of the 2D wave vector Q. The partition function 
Qn of this Coulomb fluid composed of charged particles is given by 



with 



(3) 



+ ^/3 II u(r,r')po(r)po(r)rfWr' 



(4) 
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where Po{r) = cr5(a — z) + crS{a + z) stands for the fixed charge density of 
the two walls. The three terms of /3?7[rj, Vj] correspond to direct electrostatic 
interaction between counterions, electrostatic interactions between counte- 
rions and fixed charges and between fixed charges on the walls themselves, 
respectively. 

We now proceed in the Netz [6] fashion and perform the Hubbard-Stratonovitch 
transformation of the partition function [7], where the configurational inte- 
gral over counterion positions is transformed into a functional integral over 
a fiuctuating field 0, which can be interpreted as the local fiuctuating elec- 
trostatic potential. In this way the grand-canonical partition function is 
obtained straightforwardly in the form [12] 



= 

N=0 



oo „ 

X;A^g;v = C /P[0(r)]e-^^[^W], (5) 



where A is the bare fugacity which is the exponential of the chemical potential. 
The prefactor C above is the functional determinant of the inverse interaction 
potential M~^(r, r') while the "action" of the functional integral is given by 

PH[(P] = ^pjje{r){V(P{r)fdh - ipjpo{r)(P{r)dh - \'je'^'°'''^^''^dh. (6) 

Here we have denoted the renormalized fugacity as A' = A exp(ieQg^uo(r, r)). 
The functional integral Eq. [5] can be further rescaled [12] yielding Zq — > 

C J T)[(j)(r)] e ^1^, with dimensionless H'[(j){r)], showing explicitly the de- 
pendence on the coupling parameter S. The difference between the case with 
and without dielectric discontinuities is thus only in the spatial dependence 
of the dielectric constant indicated by e(r). In this field-theoretical represen- 
tation the grand canonical partition function is now ready to be evaluated 
explicitely and analytically in the limiting cases of the weak and the strong 
coupling limits, as defined in the previous section [2,6]. 



3 Weak coupling limit 

In the limit of H ^ 1 the functional integral Eq. [5] is dominated by the 
saddle-point value of the fiuctuating potential, ^olr), set by the extremal 
value of H[(j)] [6, 7, 12]. This is defined as a solution of 

6H 



= 0. (7) 

00 
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The lowest order correction to the saddle-point is obtained by expanding 
around the extremum 0o(r) to the second order, the first order being 
zero by definition. To the lowest order in the deviation from the saddle-point 
fiuctuating potential 

= H[Mr)] + {j^^^^^ 50(r)50(r') dh dh' + 0{6ct>% 

(8) 

with 5(j){v) = 0(r) — 0o(r)- The second order term obviously describes fiuctu- 
ations around the mean-field. By making a Wick rotation in the saddle-point 
potential, (pQ — > iipo, and writing Eq. [7] explicitly, we obtain 

V2^o = -^e-^^"'?^"^. (9) 

eeo 

This is obviously nothing but the mean-field PB equation, as first shown 
in [5]. The saddle-point ansatz thus corresponds exactly to the mean-field 
approximation of classical statistical mechanics of Coulomb fiuids. In the case 
of surface distribution of external charges, as in the case of system Fig. [H the 
PB equation has to be supplemented by an appropriate boundary condition 
at z = ±a corresponding to the electroneutrality of the system. The Hessian 
of the Hamiltonian can now be evaluated explicitly and the corresponding 
grand canonical partition function for the mean-field plus the second order 
correction to the mean-field can be obtained straightforwardly as [6, 7] 

Zg = C exp(^-f3H[(j)o] - \ Trln[M-i(r, r') + A'/^e^g^ g-ZSeogVoW^s^^ _ 

(10) 

The second order fuctuational term can be shown to be equivalent to the 
first order loop expansion term [6]. As already noted, the second order fiuc- 
tuational correction to the mean-field presents a generalization of the zero- 
frequency Lifshitz interaction term [5]. In what follows we shall now evaluate 
Eq. [10] explicitly for the model system Fig. [H 

3.1 Mean- field 

The solution of the Poisson-Boltzmann equation Eq. [9] in the geometry of 
Fig. [T]is well known (see e.g. Ref. [1]) being given by iI)q{z) = In cos az 
with = X' f3eQq^ /2eeo, that can be determined from boundary conditions 
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at the two bounding surfaces z = ±a. They are of the form: atgaa = 
where /i is the Gouy-Chapman length. From the saddle-point solution 
one can obtain the corresponding grand canonical theormodynamic potential 
q = —kT In Zg- It implies a free energy defined by Legendre transform F = 
q + flN, where jl represents the chemical potential fl = kT In a^. Expressing 
all quantities in dimensionless units, so that d = a/fi, a = a/fi and p = 
p/{a'^/2eeo) the mean-field free energy is obtained as [1] 

Fo/S = 2a'^d + 2ln{a^ + 1) . (11) 

Here we have discarded all terms that do not depend on the intersurface 
separation a since we are interested solely in the interactions between the 
two charged surfaces. In the above free energy we also omitted a term that 
would go as the logarithm of C in the grand canonical partition function Eq. 
[TOl This term actually leads to an additive contribution to the free energy Eq. 
[TT] that is however exactly cancelled by the zero-frequency Lifshitz - van der 
Waals term for two dielectric half-spaces that one has to add to the total free 
energy in a consistent treatment of the interactions (for details see [5, 11]). 

The pressure corresponding to the free energy Eq. [11] and acting between 
the two surfaces is obtained simply via a derivative of the free energy with 
respect to intersurface volume, po = —l/S{dFQ/d{2a)). In dimensionless 
form it can be written as [6] 




We see that on the mean-field level the pressure is always repulsive and 
does not depend on the dielectric discontinuities in the system. This is a 
simple consequence of the fact that the mean-field solution depends only on 
the longitudinal coordinate z, whereas image effects always depend on the 
induced charge density at the interface that necessarily has also transverze, 
i.e. {x,y), dependence. On a consistent mean-field level one thus can not 
expect any image effects. 

3.2 Second order correction 

We now turn to the second order expansion term, i.e. the second term in Eq. 
[TOl We will deal with it only briefly since it was already derived before (see 
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Refs. [5, 11]). From the saddle-point solution we first of all get 



A'/3e2g2 e-P-o,Mr) = (13) 

cos^ az 

The tracelog in Eq. [TO]can be written equivalently in the following form [5,11] 

Trln[M-i(r, r') + X'(3ey e-^^°^^''('-)(53(r - r')] = — T Q dQ. (14) 

27r Jo Vq[Q) 

The discrete sum of eigenvalues of the operator + V has been replaced 
by an integral over the transverze wave- vector Q = {Qx-, Qy) with density of 
modes S'/(27r)^, due to the translational homogeneity of the system in the 
transverse directions, and Vx is the secular determinant of the fluctuating 
modes. The above formula is obtained by using the argument principle [11]. 
The index A here refers to the eigenvalue equation that can be derived in the 
form 

^ g^-A^^)/.(Q,.) = 0. (15) 



Kdz"^ cos^a^;. 
The linearly independent solutions of this equation for A = 1 inside the region 
z < \a\ are yi{z) = e^^(l + ^tgaz^ and 1/2(2) = ^^2(1 - ^tga^;). For 
the external region z > \a\ the solutions are simply exponential functions 
g±Qz_ rjj^g secular determinant is then obtained in the following form 

^i(Q)=^7(Q)(l-A2(g)e-^Q'^), (16) 



with 

9{Q) 



[(l + A)(l+aV) + 2g/i + 2g V' 
g2 + a2 



2, ,212 



^ (i + A)(i + aV) + 2A(gV-g/.) 
(i + A)(i + «v) + 2(gv + g/i) ■ 

The dielectric jump at the boundaries of the system z = ±a is characterised 
by A = (e — e') /{e + e') . We can now finally write down the second order 
correction term F2 to the free energy in dimensionless units as 

F2/S = \'E / Q\ng{Q)dQ+ Q \n{l - A^Q) e-^"^^) dQ . (18) 
'-Jo Jo 



9 



which is obviously proportional to the coupling parameter S. In order to 
extract the finite contributions to the free energy, we have to regularize the 
first integral in the above expression that yields 

(19) 



f Q\ng{Q)dQ = l-di'^{A+\n2a^)-h{&) &Tctgh{a)-—^ ^^In 

io 2 4 



where 



= VI + 2A + 2(1 + A)a2. (20) 

The second integral in Eq. [18] is convergent and can be straightforwardly 
evaluated numerically. 

Before actually proceeding numerically we first consider the limits of small 
and large spacings analytically. The second-order free energy Eq. [18] allows 
for a simple asymptotic expression valid for small d of the form 



F2/Sr.-E^^. (21) 



^ Li3(A^ 
32 a2 

Here, Lis is a polylogarithm function. This contribution comes from the 
second term of Eq. [18] and is valid for A 7^ 0. The first term of Eq. [T8] 
contributes only to the order 1/d and is therefore not inculded in the asymp- 
totic expression. The corresponding pressure is obtained from the derivative 
of the free energy with respect to the intersurface volume, and if we combine 
it with the mean-field expression we get for the total pressure in the a <C 1 
limit 

1 1 _ ^ Li3(A2^ 

d ~ 3 ^ " 32 53 

Here one needs to remember that the whole approach, meaning mean-field 
plus quadratic fiuctuations around the mean-field, only works if the fiuc- 
tuation contribution is a small perturbation to the dominating mean-field 
contribution. Thus the above formula can be assumed correct only if the 
last term is smaller than the first two. Note that only the fiuctuation term 
depends on the dielectric discontinuity. 

In the opposite limit of large a, the asymptotic form of the dimensionless 
second order free energy can be derived in the form 



P = Po+P2~^---S 00-3 • (22) 



32 d 



2d^ . „ 1 + A 



^2/^ ~ 27rMln— + 1 + 2 arctgVl + 2A 



VI + 2A 



1 + A)ln(l + A) +C(3) . (23) 
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The largest contribution goes as Ina/a^ and does not involve image effects, 
i.e. is independent of A. Image contribution can only be discerned at the 
order \jc? and its magnitude increases with increasing A. The asymptotic 
expression for the total pressure for a S> 1 is then 



P = P0+P2 ~ -TX 



32 53 L 



27rMln 



2a 



7.2 



1 + A 

I 

Vl + 2A 



arctgVl + 2A 



:i + A)ln(l + A) +C(3) 



(24) 



In this limit the fluctuation terms are consistently smaller than the mean-field 
term. Again only the fluctuation term depends on the dielectric discontinuity. 



3.3 Numerical results 




Figure 2: Rescaled pressure between the charged walls located aX z = ±a, 
as a function of the rescaled half-distance a. The mean-field approximation 
is compared with the total pressure involving the second-order correction for 
various dielectric jumps A and for coupling parameters S = 0.5 and 1. 



Numerical results for the total pressure acting between the walls, corre- 
sponding to the mean-field plus quadratic fluctuations around the mean-field, 
and given in the dimensionless form by 

p{a) =po{d) +p2{a), (25) 
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and are presented in Fig. [2] as a function of dimensionless separation a. The 
second order correction p2{a) obviously gives an attractive contribution to 
total pressure. As P2{a) becomes comparable to the leading order mean-field 
term po{d) the loop-expansion is expected to break down. As a consequence, 
this method can not be used to predict a net attraction between two charged 
plates. Since only the fluctuation term, P2(5), depends on the dielectric 
discontinuity its effects are by necessity small in this limit. For large coupling 
parameters S a different, non-perturbative theory must be used to which we 
turn next. 



4 Strong coupling limit 

In the strong coupling limit, where 2^1, the saddle-point expansion used 
before does not work [2]. Here, a straightforward perturbation expansion in 
terms of fugacity A can be used instead, corresponding formally to the virial 
expansion of the grand potential Eq. [5l For our purposes only the first two 
terms of this expansion will be used, leading to 

Zg = 4°^ + A'zg) + C(A'2). (26) 

Here the zeroth order term corresponds to the bare electrostatic interac- 
tion energy of charged walls without any counter-ions, and the first order 
term corresponds to a one-particle contribution [6]. By making an inverse 
Hubbard-Stratonovitch transformation, the zeroth order term reduces to a 
very simple form 

On the other hand, the first order term is a bit more complicated, but nev- 
ertheless leads straightforwardly to 

Z^G^ = Z^G^ y g-/3eog/«(r,R)po(r)d3r-i/3e2g2«(R,R)^3j^_ (^28) 

Note here that in both cases we are dealing with the complete electrostatic 
interaction potential, including the image term Eq. [TJ For simplicity we 
now introduce U, so that the grand partition function can be written as 
Zg = Z^il + XU), with 



U = y exp (^-^pelq^Uirn{'R,'R) - Peoq j n(r, R)po(r)d'r^ rf^R 



(29) 
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To calculate the pressure between charged walls, we now follow the same 
procedure as we did in the saddle-point approximation. The first step is to 
evaluate the grand potential q from the grand partition function Zq- Be- 
cause of the relation \U <C 1 applicable in this limit, we can expand the 
corresponding logarithm and the grand potential can be simply written as 
q = Qo — kTXU, where we have also introduced the zeroth order grand po- 
tential go = —kT In Zq\ Free energy can again be obtained by adding Gibbs 
free enthalpy flN to the grand potential. Since free energy does not explicitly 
depend on chemical potential i.e. fugacity, dF/dX = leads to a relation 
between the number of counter-ions and \U of the form = XU. The free 
energy, up to an irrelevant constant, is finally obtained in the strong coupling 
limit as 

F{a) = go (a) - NkT In U{a) . (30) 

In order to evaluate the integrals go and U, we can first integrate out the 
lateral coordinates, since our system is translationally homogeneous in the 
lateral directions. The zeroth order is then given by 

go = I dzdz'u{Q, = 0,z;z')po{z)po{z'). (31) 



Inserting the Green function u(Q, z, z') and the charge density of the walls 
Po(r) = o'6{z — a) + a6{z + a) the above expression is transformed into 

^0 = --—^ • (32) 



Here we have discarded all therms that do not depent on a since again we 
are interested only in the interactions between the two charged walls. In the 
same fashion we evaluate also the two terms involved in expression U, Eq. 
The self-image term can be expanded in series of A and then integrated 

2/ 



f d Q 

'Uim(R,R) = J j^--:^Ui,,^{Q,z;z) 



" (33) 



Aireeo ^ rfic? — z^ Svreeoo 

n odd 

Here the sum is over all odd integers n = 1,3, 5,.... The other term can be 
expressed in a much simpler form, namely 

d^r«(r,R)po(r) = -^(i±^y. (34) 
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The dimensionless free energy is finally obtained within the strong coupling 
approximation as 



2 „ „1 



F/S = 2a[l^^\ - ^ln(l - A") -21na / exp 



n odd 



(35) 



From this dimensionless free energy the dimensionless pressure is obtained 
analogously to the weak coupling case. We note here that in the case of no 
images, corresponding to a dielectrically homogeneous case and thus to A = 
0, the dimensionless pressure can be evaluated explicitly giving p = 1/a — 1, 
as shown in [2, 6]. 

It is again interesting to look at the asymptotic behaviour of pressure 
for small and large wall separations. However, one needs to exercise some 
caution here. Let us first assume that A > and is finite. In the case of 
a ^ 1 the main contribution to the integral Eq. [35] comes from small values 
of the integration variable t. Therefore, we can expand the rational function 
of t to the second order and obtain a Gaussian integral. The asymptotic 
expansion to the lowest order for a <^ 1 is then obtained as 

F/^(5,A >0) ~2afi^^ - -ln(l- A) -31n5 + C(a2), (36) 
\ 1 — Ay a 

with the corresponding pressure 

1 + A\^ S 3 



p(a.A>0)~-^— j -_,„(!- A) + - + 0(a). (37) 

The interaction pressure obviously depends strongly on the value of the di- 
electric discontinuity. As A — > one has to carefully pick the terms in Eq. 
that remain finite, ending up with 



F/S{d,A — > 0) = 2{a- In 5) + (8a + =) A + ^(A^). (38) 
with the pressure 

A^0) = ^- l + (^-4)A + OiA'). (39) 

Clearly the pressure for A = reduces top = 1/a — 1, which is consistent with 
the result obtained with no images [2,6]. The above two limits do not imply 
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any discontinuity for pressure at A = 0. The change from the asymptotics 
Eq. [37] to [39] is smooth and continuous in A. 

In the opposite asymptotic hmit of large a, the exponential function in Eq. 
[35] is approximately 1 and the free energy can be written straightforwardly 
as ^ 

F/S r^2d(\^^ - ^ln(l - A^) -21na + 0(a-2). (40) 
The corresponding pressure is then 

Again the dependence on the dielectric discontinuity is quite strong. The 
expression involving no dielectric discontinuities is here obtained simply by 
setting A = directly in Eq. 14T1 

Apart from the interaction pressure we also evaluate the counterion den- 
sity in order to quantify the image effects on the distribution of mobile charges 
between the two surfaces. The counter-ion density can be evaluated follow- 
ing [6] by using the relation N = XU and inserting the expression for U. The 
integrand can then be recognised as the counter-ion density of the form 



p{z) = A{a) exp 



n odd 



(42) 



The normalization factor A{d) can be obtained by the volume integration of 
density and the electroneutrality ansatz. 

At the end we note here, that we have to be aware of the range of va- 
lidity of the strong coupling limit (S ^ 1) which is formally given by the 
Rouzina-Bloomfield criterion [13] H > as formulated by Netz [6]. At a 
given coupling parameter the strong coupling limit thus applies only to suf- 
ficiently small intersurface separation. Conversely at a given separation the 
strong coupling limit applies only to sufficiently high values of the coupling 
parameter. 



4.1 Numerical results 

We now present numerical results for the dimensionless interaction pressure 
in the strong coupling limit. It depends only on a and S. In the general 
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Figure 3: Rescaled interaction pressure between the charged walls as a func- 
tion of the rescaled half-distance a for coupling parameter S = 10 and various 
values of the dielectric jump defined by A. 

case of nonvanishing A the pressure has to be evaluated numerically and is 
shown in Fig. [3l It is obvious from figure Fig. [3] that the rescaled pressure 
depends strongly on the value of the dielectric discontinuity A. For larger A 
the attraction at larger separations as well as repusion at smaller separations 
are larger. The attractive pressure obviously diverges as A approaches its 
limiting value of 1, valid for ideally polarizable surfaces. 

We note that just as in the case of no dielectric inhomogeneities [12] 
there exists an intermediate value of the intersurface separation a* at which 
the pressure equals zero. This separation defines the bound state of the 
system and has been studied extensively [12]. In the case of dielectric inho- 
mogeneities the position of the bound state depends drastically on the value 
of A, see Fig. HI For A = the bound state is at a* = 1 [12]. This value in- 
creases with increasing A and reaches a maximum, then decreasing down to 
zero as A approaches 1. This non-monotonic dependence of the bound-state 
on the dielectric discontinuity is quite interesting. In the case of an ideally 
polarizable interface in the strong coupling limit we would thus have com- 
plete collapse with d* = 0. Here the correlation effect obviously overwhelms 
the entropy of the conterions, collapsing the system completely. We just note 
here that for A < the strong coupling pressure is purely attractive and the 
system would collapse completely for any A < 0. 

While the dielectrically homogeneous system has been studied extensively 
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Figure 4: Rescaled half-separation between the walls of bound state as func- 
tion of dielectric jump A for various values of coupling parameter S. 

and the analytic calculations can be exhaustively compared with existant 
computer simulations this is not the case for dielectrically inhomogeneous 
systems. The presence of image effects and the underlying image potential 
restricts the computer time considerably and the corresponding simulation 
studies are scarce. The only available computer simulation data are the 
Monte Carlo simulations by Naji et al. [2]. Unfortunately even these simula- 
tions can not be compared directly with our results. The main deficiency of 
these simulations is the approximate description of the image effects, since 
only the first-order images are taken into account. Focusing on the bound- 
state wall separations for A = 0.95 (corresponding to water-hydrocarbon 
interface), the simulation results are an order of magnitude larger than our 
numerical results. The major part of discrepancy is surely caused by the 
first order approximation in the simulations. Namely, if we expand our ex- 
pression for the free energy Eq. [35] to the first order in A, thus making a 
similar approximation as in the simulations, our results are only about 30 % 
off, as shown by dashed line in Fig. [51 From this we can conclude that image 
contribution to the interaction is of crucial importance. 

At the end let us also analyze results for counter-ion density as shown 
in Fig. [6l We note that in general the image effects change the distribution 
of counterions. For repulsive image interactions the counterions tend to be 
depleted in the vicinity of the charged surfaces. We note that this depletion 
appears to be weak for small dielectric jumps A and becomes significant for 
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Figure 5: Bound state in terms of separation a* as a function of H at di- 
electric jump A = 0.95 at both walls. Dashed line represents the linear 
approximation and circles correspond to Monte-Carlo simulations for first- 
order images [2]. 
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Figure 6: Density profile between two charged walls separated by distance 
2a = 2 for various dielectric jumps A at coupling parameters S = 1 and 
H = 10. Here the volume integral is normalized to 1. 



larger A. This of course makes sense, since higher A means larger charge of 
images on the other side of the wall that causes higher electrostatic repul- 
sion, so counter-ions tend to concentrate toward the center. This effect also 
increases with increasing coupling constant S. Obviously, the equation Eq. 
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H2] is reduced to a constant when A vanishes, since in that case there is no 
image repulsion. 

5 Conclusions 

In this paper we have analysed the effect of dielectric discontinuities on the 
interaction pressure between two charged surfaces with intervening counteri- 
ons in planparallel geometry. In this analysis we have followed the approach 
of Netz and coworkers by explicitly considering the limits of weak and strong 
coupling. This limits refer to the value of the electrostatic coupling param- 
eter S, which is assumed small in the first limit and large in the second, 
respectively. 

Consider first the weak coupling limit with additive contributions from 
the mean-field and harmonic fluctuations around the mean-field. We note 
that the mean-field term in the interaction pressure is not affected by image 
interactions. This is easy to understand: in planparallel geometry the mean 
field depends only on the transverse coordinate, while image interactions 
by necessity imply also dependencies in the transverse directions. Thus by 
construction the mean field solution in this geometry does not depend on the 
presence of electrostatic image interactions. One could add those in by hand, 
as was done recently by Onuki [14], but such an approach does not sensu 
stricto correspond to a mean-field analysis. On the weak coupling level the 
image interactions thus only transpire in the second-order fluctuation term 
around the mean-field and have the form of a generalized zero-frequency 
Lifshitz term. Since this term in the weak coupling limit is by necessity 
small, so is in general the image effect. 

The opposite is true for the strong coupling limit which is effectively 
a single particle description since it corresponds to the second order virial 
expansion. As such it is quite susceptible to dielectric inhomogeneities and 
their effect on the interaction pressure is substantial. The image effects boost 
the value of the interaction pressure in the repulsive - small separations - as 
well as attractive -larger separations - domain. For the infinitely polarizable 
case in fact the pressure can reach infinite values. Also the bound state, 
defined as the separation at which the interaction pressure is zero, depends 
drastically and non-monotonically on the value of the dielectric discontinuity. 
For large differences in the value of the dielectric constant the system shows 
a tightly bound state. 
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Our computation has an important practical aside. As we have shown, the 
first order image correction is not enough to describe accurately the dielectric 
discontinuity effects on the strong coupling level. Even simulations that take 
into account only the first order images are not particularly accurate and 
miss most of the image effect in the interaction pressure. Our calculation 
is thus particularly valuable in the region of parameter space where other 
methods obviously fail quite seriously. 

The validity of our approach is set by the constraints of the weak-strong 
coupling formalism as well as its coupling to the image effects. This is par- 
ticularly clear in the latter. The approach followed in this contribution only 
makes sense if effective lateral separation between counterions is larger than 
the separation between images in the longitudinal direction. If this is not the 
case the strong coupling limit as concieved of above could not be constructed. 
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